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ABSTRACT 

A  Nusselt-type  analysis  was  performed  for  laminar  film  con- 
densation on  the  inside  of  a  truncated  rotating  cone.   Studies 
of  this  model  revealed  that  the  heat  transfer  capacity  of  the 
rotating  thermosyphon  was  controlled  by  three  thermal  resis- 
tances; the  condensate  film  thickness,  condenser  wall  resistance, 
and  the  external  heat  transfer  coefficient.   A  parametric  study 
of  the  solution  was  performed  for  various  values  of  these 
resistances  and  methods  for  reducing  them  were  investigated. 

A  stainless  steel  thermosyphon  was  tested  using  various  work- 
ing fluids.   A  copper  thermosyphon  was  constructed  and  tested 
using  water.   Different  power    levels  and  speeds  of  rotation 
were  investigated.   Efforts  were  made  to  increase  the  external 
heat  transfer  coefficient  by  regulation  of  the  cooling  spray. 

Longitudinal  grooving  was  studied  as  a  means  of  reducing 
the  heat  transfer  resistance  of  the  internal  condenser  surface. 
An  analytical  model  was  developed  for  a  triangular  groove  pro- 
file and  a  parametric  analysis  performed  to  show  effects  on 
heat  transfer  capabilities. 
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I.   INTRODUCTION 


A.   THE  ROTATING  TWO-PHASE  THERMOSYPHON 

The  rotating  two-phase  closed  thermosyphon  is  a  device 
capable  of  transferring  large  quantities  of  heat  from  one 
place  to  another  at  nearly  isothermal  conditions.   The 
three  main  parts  of  the  thermosyphon  are  the  evaporator, 
condenser,  and  the  working  fluid  (Figure  1). 


MEAT    IN 


HEAT      OOT 


o 


Fig.  1   Schematic  Drawing  of  Typical  Rotating  Two-Phase 
Thermosyphon 


Working  fluid  is  evaporated  in  the  evaporator.   A  small 
pressure  gradient  forces  the  vapor  to  flow  towards  the  con- 
denser where  the  latent  heat  of  the  vapor  is  removed  upon 
condensation.   The  condenser  wall  has  a  small  degree  of  taper 
opening  towards  the  evaporator.   A  high  centrifugal  force 
created  by  rotation  of  the  thermosyphon  assembly  about  its 
longitudinal  axis  then  drives  the  condensate  back  to  the 
evaporator.   It  is  this  external  force  field  for  returning 
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the  condensate  that  distinguishes  the  thermosyphon  from  a 
heat  pipe.   A  heat  pipe  employs  wicking  material  to  provide 
the  necessary  driving  force  for  the  condensed  working  fluid. 
The  thermosyphon  may  be  stationary  and  depend  on  earth 
gravity  to  provide  this  driving  force  or  it  may  be  in  motion 
and  rely  upon  accelerations  (centrifugal  or  otherwise)  for 
the  force.   Because  there  is  no  wick  in  the  thermosyphon 
and  therefore,  no  large  wick  flow  resistance,  the  thermosyphon 
can  transfer  large  quantities  of  heat  [1,  2]. 

B .   BACKGROUND 

Ballback  initially  studied  the  operational  principles  and 
performance  characteristics  of  the  rotating  thermosyphon  [3]. 
His  analytic  development  was  a  Nusselt-type  analysis  for  film 
condensation  on  the  inside  surface  of  a  rotating  truncated 
cone.   He  obtained  a  closed  form  solution  for  condensate 
film  thickness,   &   ,  by  incorporation  of  the  following 
assumptions  in  his  analysis:   (1)   the  rate  of  change  of  film 
thickness  along  the  x-axis  of  the  condenser  (Figure  2)  was 
much  less  than  the  tangent  of  the  half  cone  angle,  2-^«tan  4>  , 
(2)  there  were  no  shear  stresses  between  the  vapor  and  con- 
densate, and  (3)  there  was  no  thermal  resistance  in  the 
condenser  wall. 


.  J 

Figure  2   Coordinate  System  and  Geometry  for  Condensate 
Film  Analysis 


Daley  [4]  performed  a  subsequent  analysis  which  included 
the  slope  of  the  film  thickness  and  thermal  resistances  due 
to  the  consenser  wall  and  outside  cooling  mechanism.   A 
second  order  differential  equation  for  the  condensate  film 
thickness  resulted.   This  equation's  non-linearity  proved  to 
be  very  sensitive  to  the  initial  value  of  film  thickness  and 
consequently,  no  solutions  were  obtained  for  half-cone  angles 
greater  than  zero.   Newton's  study  [5]  improved  the  analysis 
by  including  interfacial  shear  stress  between  vapor  and 
condensate  and  accounting  for  change  of  vapor  pressure  along 
the  condenser  length.   The  problem  was  reduced  to  two  simul- 
taneous first  order  ordinary  differential  equations  having  as 
unknowns  the  condensate  film  thickness,   &(x),  and  the  vapor 
velocity,  v (x) .   A  Runge-Kutta  numerical  integration  scheme 
was  used  to  solve  these  equations  on  an  IBM  digital  computer. 
The  integration  was  started  by  assuming  initial  conditions, 
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i   =    i,     ,  v   =  0,  and  p.   =  taM  at  x  =  °'   Results  were 
1  d  x 


obtained  for  half-cone  angles  of  0°,  0.1°,  0.2°,  and  0.3°. 

This  solution  showed  £L±  «  tan  <t>  and  allowed  simplification 

dx 

of  the  governing  equation  for  &  (x) .   The  approximate  solu- 
tion deviated  from  the  numerical  solution  by  a  conservative 
2.5%  and  agreement  improved  as  the  half  cone  angle  was 
increased.   For  half  cone  angles  greater  than  0.3°,  the 
approximation  gave  solutions  that  were  independent  of 
reasonably  small  values  of  initial  film  thickness,  d\_   . 
Woodard  [6]  followed  this  development  with  a  parametric  study 
of  the  thermosyphon  operating  characteristics  using  the  approx- 
imate solution  of  Newton. 

Daley's  study  concluded  the  design  and  assembly  of  the 
thermosyphon.   Newton  operated  the  apparatus  and  was  able  to 
make  the  first  comparisons  between  analytical  and  experimental 
results.   Experimental  data  was  2  0%  higher  than  maximum 
theory  due  to  uncertainty  of  parameter  measurements  and  exis- 
tence of  dropwise  vice  film  condensation.   Woodard  was  able 
to  improve  the  instrumentation  and  implement  thermosyphon 
cleaning  procedures  that  promoted  film  condensation.   He  then 
obtained  additional  data  for  comparison  with  theory  and  also 
found  that  the  condenser  wall  temperature  was  not  constant 
along  its  length. 

C.   THESIS  OBJECTIVES 

The  objectives  of  this  study  were  to:   (1)   investigate 
the  previous  analytical  solution  in  an  attempt  to  improve 
thermosyphon  performance,  (2)  improve  instrumentation  in 
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order  to  achieve  more  reliable  experimental  data,  (3)  obtain 
experimental  results  from  different  test  sections  for  compar- 
ison with  theory,  (4)  develop  a  theoretical  model  for  a 
grooved  condensing  surface  aimed  at  increasing  the  rate  of 
heat  transfer,  and  (5)  conduct  a  parametric  study  of  the 
grooved  versus  smooth  condensing  surfaces. 
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II.   THEORETICAL  PROGRAM 

A.   REVIEW 

To  improve  the  Nusselt-type  analysis  performed  by  Ballback 
and  Daley,  Newton  included  the  effects  of  interfacial  shear 
stress  between  liquid  and  vapor,  and  vapor  pressure  change 
along  the  condenser  axis.   Momentum,  continuity,  and  energy 
equations  were  developed  for  an  infinite  simal  condensate 
element.   Momentum  and  continuity  equations  were  also  devel- 
oped for  one-dimensional,  turbulent,  and  isothermal  vapor 
flow.   Reduction  of  these  equations  in  Newton's  study  [5] 
resulted  in  two  simultaneous  first  order  ordinary  differential 
equations 


dv    =     2  (Ro+-xs.n4>KT5-Tj     _  zv  r  ^  ,    _  Cos  x  i£  1 


i) 


r  r  p>     Pvi  r^k?     o  d  dv  i 

ii    =    tan  <j>  -    pv^LR^Tj^    PvvL-fe    lUil] 
dx  P,cos<t>    [pf^R  +.   l^    ] 

P,    =      (Ro+^m«f))|     -     COS     (j>      A     S" 

Pt  -     (R0+-X  sin  4>;A.    ~    cos  dj>     -£- 

with  the  following  initial  conditions: 

£(o)  =    Si 
v(o)  =  o 


(2) 
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A  Runge  Kutta  numerical  integration  scheme  was  used  to  solve 
these  equations  on  an  IBM  digital  computer.   Solutions  were 
highly  dependent  upon  the  choice  of  initial  film  thickness 
and  because  of  instability  problems,  were  only  obtainable  for 
half  cone  angles  of  0°,  0.1°,  0.2°,  and  0.3°.   The  results 
indicated  that,  except  near  the  ends  of  the  condenser, 

r\  i. 

— -  «  tan  <b    .   This  approximation  and  an  order  of  magnitude 
dx 

analysis  applied  to  equations  (1)  and  (2)  reduced  them  to  a 
first  order  differential  equation  in  v (x)  and  a  cubic 
equation  for  &  (x) 

dv  _   2(Ro+X  Sin  <j>)  (Ts-Tcq)  Zv  Sin  4> 


•£rsm  +  [ff  "lR  +^]-  f.v'f  [R-|]-  fvv^  =o  (4) 


Equation  (3)  with  the  use  of  equation  (4)  was  numerically 
integrated  over  the  length  of  the  condenser  using  initial 
values  of  zero  film  thickness  and  zero  velocity.   Having 
found  vapor  velocity  at  a  particular  step,  a  subroutine  was 
used  to  solve  equation  (4) .   This  resulting  value  of  film 
thickness  was  used  to  find   v(x)  for  the  next  step  of  the 
integration.   The  total  heat  transfer  rate  was  then  obtained 
by  summing  the  heat  contribution  of  each  incremented  step  over 
the  condenser  length 

2.1T  (R0  *  X  sin  4>)  (Ts  -To)  A*  ,_. 

1«  +  .1-  4-  7- 
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Results  from  equations  (3)  and  (4)  were  2.5%  lower  than  the 
more  exact  solution  (from  equations  (1)  and  (2))  for  half  cone 
angles  of  0.3°  and  less.   For  larger  half  cone  angles  (up  to 
3.0°)  the  approximation  gave  results  that  were  independent  of 
the  small  initial  film  thickness.   Thus,  equations  (3),  (4), 
and  (5)  were  used  to  obtain  theoretical  results  for  condensate 
film  thickness  and  heat  transfer  from  the  thermosyphon. 

B.   MODIFICATIONS  TO  THERMOSYPHON  THEORETICAL  RESULTS 

Woodard ' s  experimental  data  indicated  that  modifications 
to  the  computerized  solution  of  equations  (3)  and  (4)  were 
necessary.   The  outside  condenser  wall  temperature  had  orig- 
inally been  taken  as  constant  over  its  entire  length.   This 
value  was  the  average  of  three  or  four  surface  thermocouples 
which  were  placed  at  different  locations  along  the  condenser. 
The  thermocouple  closest  to  the  evaporator  usually  indicated 
a  higher  temperature  than  the  remaining  thermocouples.   Re- 
location and  addition  of  thermocouples  in  subsequent  runs 
established  the  existence  of  a  steep  temperature  gradient  on 
the  condenser  exterior  near  the  evaporator.   Equations  (3), 

(4),  and  (5)  and  the  computer  solution  were  then  modified  to 
account  for  this  variation  in  wall  temperature.   Equations 

(6),  (7),  and  (8)  reflect  these  changes. 


dv     _     2.(rc+  x  sin  4>)  (Ts-TwCx))             2vsm$ 
_    _     _ (6) 

dx 


*[S?  + 1!*1*, 
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i?R  s.n*  [ft:  u>*R    *■  *£**]  -  PvV*f  [R  f  ]-  frV?^  =  0       (7) 
3  LX  R      J       ^  ^  v        2.  (V 


_   2LTT  (R0  +  X  5m  fr)  (Ts-Tw(x))  ax 


Since  varying  temperatures  were  encountered  not  only  from 
one  data  point  to  another  but  also  along  the  length  of  the 
condenser  as  described  above,  the  change  in  fluid,  vapor,  and 
wall  properties  with  temperature  had  to  be  accounted  for.   A 
least-squares  polynomial  was  generated  for  the  various  prop- 
erties over  the  temperature  ranges  involved.   These  polynomials 
were  placed  in  the  computer  solution  such  that  new  properties 
were  computed  at  each  incremental  length. 

The  last  modifications  made  to  the  computer  solution  in- 
volved the  friction  factor  used  to  account  for  shear  forces 
at  the  vapor-condensate  interface.   A  Reynolds  number  for  the 
vapor  flow  was  calculated  at  each  incremental  condenser  length. 
Engineering  judgment  indicated  that  turbulence  in  the  vapor 
flow  would  occur  at  approximately  the  same  critical  Reynolds 
number  as  pipe  flow.   Accordingly,  the  computer  solution  was 
changed  to  use  laminar  or  turbulent  flow  friction  factors 
when  the  vapor  flow  Reynolds  number  was  above  or  below  a 
critical  value  of  2300.   This  effect  was  found  to  be  small 
as  evidenced  by  results  obtained  using  a  friction  factor  ten 
times  that  normally  encountered. 
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C.   PARAMETRIC  STUDY  OF  SMOOTH  CONDENSER 

Woodard ' s  thesis  [6]  included  a  parametric  study  for  the 
thermosyphon  with  a  smooth  stainless  steel  interior  surface. 
The  variables  considered  were  condenser  geometry,  rotational 
speed,  working  fluids  and  the  exterior  heat  transfer  co- 
efficient.  In  an  attempt  to  improve  the  heat  transfer 
capabilities  of  the  thermosyphon,  it  became  necessary  to 
determine  those  factors  which  limited  performance.   Represen- 
tative values  of  condensate  thickness  and  average  external 
heat  transfer  coefficient  taken  from  Woodard ' s  experimental 
results  gave  the  following  thermal  resistances: 


THERMAL        RES\STANCE" 

REPRESENTATIVE        VALUES 

INSIDE 

i 

i   -Z.    .00035     (-ft) 

kf-   .2,5      (Btu/kr- ft -°F) 

-r   ~    .001 

WALL 

t 

t    ~    .OOSZ       (ft) 

kw~  °i.o  (Btu/V-ft-T) 

t 

—  z-  .oook 

OOTStDE 

HEKT     TRANSFER 

COEFFICIENT 

i 
h 

12.00     <     h„     <      IbOO 
(B-tu/U.r-ftl-0F) 

-r  X   .000^ 

Table  1   Thermosyphon  Thermal  Resistances  from  Data  of 
Woodard  [6] 


Table  1  clearly  indicates  -  the  relative  importance  of  each 
thermal  resistance  and  shows  those  factors  having  strongest 
influence  on  overall  heat  transfer  capabilities.   The  external 

heat  transfer  coefficient,  h  ,  was  found  to  vary  from  1200  to 

2 
1600  Btu/hr-ft  -°F  and  thus  rival  the  condensate  film  thick- 
ness (about  .00035  ft)  for  dominance. 
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1 .  Wall  Resistance 

Although  not  the  controlling  factor,  wall  resistance 
is  the  easiest  to  reduce.   A  copper  test  section,  for  example, 
would  reduce  this  resistance  by  a  factor  of  25.   Of  note, 
however,  is  the  fact  that  copper  represents  the  material 
possessing  the  highest  thermal  conductivity  other  than  silver 
and  is  therefore  a  limiting  value.   Reduction  of  wall  thick- 
ness becomes  impractical  beyond  certain  strength  limitations 
or  when  internal  grooving  of  the  condenser  is  considered. 

2 .  Internal  Resistance 

Of  ordinary  fluids,  water  possesses  the  highest  latent 
heat  of  vaporization  and  thermal  conductivity.  This  choice  of 
working  fluid  limits  reduction  of  the  film  resistance  to  some 
means  of  reducing  the  condensate  film  thickness.   Longitudinal 
grooving  of  the  interior  condenser  surface  as  discussed  later 
would  be  an  effective  way  of  doing  this. 

3 .  External  Resistance 

Lowering  of  the  external  thermal  resistance  involves 
an  increase  in  the  external  surface  area  and/or  an  increase 
in  the  heat  transfer  coefficient.   These  appear  to  be  the  most 
difficult  modifications  to  implement.   The  addition  of  small 
annular  fins  to  the  exterior  surface  does  increase  the  cooling 
surface  area  but  also  complicates  the  temperature  distribution 
in  the  condenser  wall  (an  effect  that  is  likewise  experienced 
with  internal  grooving) .   Thermocouple  location  and  the 
positioning  of  coolant  spray  nozzles  becomes  more  difficult. 
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As  a  means  of  comparison,  external  annular  fins  of 
copper  having  a  width  of  1/16  inch,  separation  of  1/16  inch 
and  protruding  from  the  surface  by  3/16  of  an  inch  were 
chosen.   Using  Holman  [7],  a  fin  efficiency  of  about  0.80  and 
overall  efficiency  of  0.82  were  obtained.   The  overall 
efficiency  and  ratio  of  areas  (total  finned  and  surface  area 
to  the  surface  area  without  fins)  were  then  used  to  obtain  an 
effective  value  of  heat  transfer  coefficient  to  use  in  the 
computerized  model  according  to  the  formula  [8] 

heff  =   ho"^  -£■  (9) 

where  hQ  =  1500  Btu/hr-f t2-°F 

n 

This  effective  coefficient  was  about  5000  Btu/hr-f t  -°F ,  or 
over  three  times  the  experimentally  obtained  value  for  exist- 
ing condensers.   It  was  then  used  in  the  analytical  model  to 
represent  the  thermosyphon  improvement  for  addition  of 
external  fins. 

Other  methods  for  improving  heat  transfer  coefficient 
include  changes  in  coolant  spray  nozzle  size,  location, 
quantity,  and  the  coolant  flow  rate.   This  would  affect  the 
droplet  size,  distribution,  and  velocity.   Lack  of  a  suitable 
increase  in  the  heat  transfer  coefficient  by  any  of  the  above 
methods  would  then  necessitate  a  thorough  study  aimed  at 
determining  the  maximum  possible  coefficient  and  then  design 
of  a  cooling  system  that  could  approach  this  value.   Such  a 
possibility  might  be  construction  of  a  coolant  jacket  that 
would  rotate  synchronously  with  the  thermosyphon.   The  coolant 
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would  have  solid  body  rotation  with  the  thermosyphon  while 
axial  flow  would  provide  the  desired  convection  cooling. 
4 .   Analysis 

Figures   3  ,   4  ,  and   5   depict  the  results  of 
changing  the  three  possible  thermal  resistances.   All  graphs 
are  plotted  with  saturation  temperature  of  the  working  fluid 
on  the  abscissa  and  the  heat  transfer  rate  of  the  thermosyphon 
on  the  ordinate.   With  the  parameters  h,  <f>  ,  and  RPM  held 
constant^  such  a  graph  represents  thermosyphon  reponse  as  a 
function  of  fluid  properties.   Variance  of  an  individual  para- 
meter then  produces  a  family  of  curves  for  each  figure  and 
thereby  illustrates  the  effect  of  that  parameter  on  the 
response  of  the  thermosyphon.   For  the  study,  values  of  ex- 
ternal heat  transfer  coefficient  used  were  1000,  1500,  and 
5000  Btu/hr-ft2-°F.   Thermosyphon  rotational  speeds  of  1000 
and  3000  RPM  were  used.   Two  condenser  materials,  copper  and 
stainless  steel,  were  compared  as  well  as  an  external  fin 
arrangement  that  would  increase  cooling  surface  area  by  4 
times.   The  condenser  half  cone  angle  of  1°  and  working  fluid 
of  water  were  unchanged  throughout  the  study. 

D.   THEORY  FOR  INTERNALLY  GROOVED  CONDENSER 

Longitudinal  grooving  of  the  inside  condenser  surface  is 
one  method  of  lowering  the  film  thermal  resistance.   Gregorig 
[9]  has  shown  that  longitudinally  grooved  vertical  condensing 
surfaces  will  have  a  net  effective  thermal  resistance  less 
than  that  of  a  smooth  condensing  surface.   He  showed  that 
this  effect  was  caused  by  a  pressure  gradient  in  the  condensate 
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due  to  the  surface  tension  forces  of  the  liquid  and  the  local 
radius  of  curvature.   For  the  geometry  of  Figure  6,  the 
pressure  at  point  1  and  2  are  computed  to  be 


P4  -  K    *- 


p  =  P 


Z  6 


Z   £ 
R, 


Figure  6.   Grooved  Condenser  Wall 

where  the  change  in  sign  results  from  the  change  in  sign  of 
the  curvature.   Assuming  vapor  pressure,  Pv  ,  to  be  constant, 
P   is  thus  greater  than  P„  and  a  driving  force  sends  conden- 
sate to  the  bottom  of  the  groove.   Although  heat  transfer  is 
poorer  in  the  trough,  very  high  heat  transfer  occurs  near  the' 
crest  resulting  in'  an  overall  increase.   When  applied  to  the 
rotating  thermosyphon ,  the  impressed  rotational  gravity  field 
would  overshadow  the  surface  tension  induced  force.   Never- 
theless, results  would  be  similar  with  the  thin  film  occurring 
at  the  crest  of  the  grooves.   An  experimental  study  by 
Lustenader,  Richter,  and  Neugebauer  [10]  has  confirmed  the 
existence  of  higher  heat  transfer  rates  for  thin  film  conden- 
sation on  a  fluted  surface. 
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For  simplicity  of  analysis,  a  groove  profile  having  a 
triangular  cross  section  was  chosen.  An  enlarged  view  of 
the  triangular  fin  geometry  is  depicted  in  Figure  7  along 
with  appropriate  dimensions.  The  x-axis  is  normal  to  the 
y-z  plane  and  is  coincidental  with  the  fin  tip.  The  radius, 
r  ,  for  a  cross  section  is  taken  arbitrarily  as  halfway 
between  tip  and  bottom  of  the  groove  and  is  represented  by: 

r,  =  Rz  +  x  sm  <(>  (10) 

This  choice  of  radius  makes  the  centrifugal  gravity  force  con- 
stant for  the  entire  depth  of  groove  at  each  x  location.   Con- 
densate from  the  triangular  fins  collects  in  the  trough  of  the 
groove  which  then  acts  as  the  primary  path  for  return  of  the 
fluid  to  the  evaporator.   A  greater  film  thickness  in  the 
trough  allows  very  little  heat  transfer  when  compared  to  the 
thinner  film  existing  on  the  exposed  fins. 
1 .   Assumptions 

The  theoretical  model  for  internally  grooved  condenser 
was  developed  based  on  the  following  assumptions  [3] : 

a.  Film  condensation,  not  dropwise  condensation, 
occurs  in  the  condenser  section. 

b.  The  condensate  film  undergoes  laminar  flow. 

c.  Subcooling  of  the  condensate  may  be  neglected. 

d.  Momentum  changes  through  the  condensate  are  small 
and  thus  there  is  essentially  a  static  balance  of 
forces. 
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2  =  0.762  inches 

t  =  0.0375  inches 

b  =  0.025  inches 

t*  =  O.050  inches 


Figure  7   Geometry  and  Coordinate  System  for 
Triangular  Grooves 
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e.  There  is  no  interfacial  shear  between  vapor  and 
condensate . 

f.  There  is  a  linear  distribution  of  temperature 
through  the  condensate  film. 

g.  Centrifugal  force  is  much  greater  than  earth 
gravity.   The  latter  may  be  neglected. 

h.   The  condensate  film  thickness  is  much  less  than 

the  radius  of  curvature  of  the  condenser  wall. 
i.   The  fluid  density  is  much  greater  than  vapor 

density, 
j .   The  thermosyphon  is  operating  at  a  steady  state 
condition. 
2.   Mass  Flow  in  Z-Direction  (per  unit  length  in  x-direction) 
In  a  manner  similar  to  Ballback's  force  balance  on  an 
inf initessimal  fluid  element,  the  momentum  equation  in  the 
z-direction  may  be  written 

°      *     -        dP       -     of  u>'L  r,    cos  <f>    cos  o<  (11) 

^  y  <^2. 

where  — H  may  be  neglected.   Integrating  equation  (11)  from 
dz 

Te  to  0  and  y  to  &   gives 

T  =  Px  WV,  cos  <b  cos  oC  (o-y) 
z     C*   ''      T        ^    '  (12) 

d  w 


T3    -      M_ 

Y 


*      ~*       * 


Integration  of  equation  (12)  from  0  to  w  and  0  to  y  gives  an 
expression  for  the  velocity  of  the  fluid 


w 


_      ft  u?1  r,     cos    fr  cos   *       (*V-4^) 


»< 


(13) 
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which  may  be  used  to  find  an  average  velocity  using  the 
expression 

-  =  _L  (  w  J.  _    f^'tX*  *  sm  4,]  cos  <).  COS*  8    (14) 

The  mass  flow  rate  at  any  z  location  on  the  fin  may  be 
written  as 

ft-  =   A  S 


'1 


f*  *  W  (15) 


Substitution  of  equation  (14)  into  equation  (15)  then  gives 
the  mass  flow  rate  as  a  function  of  film  thickness 

PrLo}L  CRj.  +  X  Sin  <b)  COS  d)  COS  crt.     S 

nt=-e 5^ (i6) 


3 .   Mass  Flov7  in  X-Direction 

An  assumption  is  made  that  the  mass  flow  rate  in  the 
x-direction  along  the  fins  is  negligible  compared  to  the 
axial  flow  in  the  troughs.   The  trough  mass  flow  at  any  x 
location  is  therefore  the  total  mass  flow  up  to  that  point. 
Using  the  x-direction  momentum  equation  for  an  inf initessimal 
fluid  element  in  the  trough 

-  -  fro>  r,  s.n<{>  (17) 

d  y        v 

and  proceeding  in  a  manner  analogous  to  that  for  the  z 
direction,  an  expression  for  average  film  velocity  in  the 
trough  may  be  obtained 

Q  -  i_  \  udy  -   pf  co2(Rr-h  x  sin  4Q  s.h  ft  6  {1Q) 
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The  total  mass  rate  of  condensate  flow  in  the  longitudinal 
direction  in  terms  of  trough  film  thickness,  b    ,  and  x 
location  becomes 


<\ot*,     =       f*   3      A^ooh 


(19) 


— (*  I      +    i     cos*  sin  <0 


TOTM. 


^-F 


where  A.     ,  =  cross  sectional  area  of  trough 
trough 

4 .   Energy  Equation  for  Fin  Condensate  (per  unit  length 
in  the  x-direction) . 

An  energy  balance  for  a  differential  element  of  the 

fin  condensate  flowing  in  the  z-direction  yields  the  necessary 

relation  between  film  thickness  h     and  distance  along  the 

fin,  z. 

dQ  -LdM  =   MT*-TWAI,U)  dm        (20) 

1     *s    i  ^ 

The  interface  wall  temperature,  T   n . ,  is  obviously  a  function 

wall 

of  distance  z.  Since  the  film  thickness  at  the  apex  of  the 
fin  (point  A  in  Figure  7)  is  zero,  it  must  be  assumed  that 
at  the  apex  the  temperature  is  the  same  as  vapor  saturation 
temperature.  It  may  also  be  deduced  that  lines  A-F  and  C-D 
represent  insulated  boundaries  because  of  profile  symmetry. 
An  adequate  estimate  of  the  interior  wall  temperature  at 
point  C  may  be  found  by  equating  heat  fluxes  for  an  infini- 
te simal  area  at  this  point 
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—   -   _ zJ—  (2\\ 


Thus,  the  wall  temperature  varies  somehow  from  Ts  to  T   along 
the  wall  from  point  A  to  C .   Observing  that  thin  films  occur 
along  the  side  of  the  fin  when  compared  to  the  film  thickness 
in  the  trough,  it  may  be  deduced  that  most  of  the  heat  flow 
occurs  through  the  region  of  the  fin.   Also,  when  considering 
one-dimensional  flow  through  section  ABEF ,  this  heat  flow  may- 
be so  great  that  its  only  limitation  is  the  amount  that  can 
be  removed  externally  from  FE .   With  this  limitation  and  the 
lower  temperatures  existing  in  section  BCDE,  it  becomes 
apparent  that  a  two-dimensional  heat  flow  problem  exists. 
Such  two-dimensional  heat  flow  in  the  condenser  wall  results 
in  a  wall-film  interface  temperature  that  may  not  be  a  linear 
function  of  z.   Until  a  detailed  two-dimensional  analysis  is 
performed,  a  linear  relationship  is  assumed  and  the  tempera- 
ture difference  in  equation  (20)  becomes 


(22) 
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Equation  (16)  is  used  to  find  the  differential,  ^Ht,  which 
with  equation  (22)  is  substituted  into  equation  (20) .   Equa- 
tion (20) is  then  integrated  from  0  to  &      and  0  to  Z   to 
obtain  an  expression  for  film  thickness,  & 


COS  o<  qx>5    <t> 


1* 


(23). 


5 .   Energy  Equation  for  Trough  Condensate 

Even  though  the  condensate  in  the  trough  has  a  larger 
thickness,  there  is  a  small  amount  of  vapor  condensation  on 
its  surface.   The  energy  balance  gives  an  expression  for  the 
mass  rate  entering  the  trough  (per  unit  length  in  x  direction) 
in  this  manner: 


dx         (24) 


where   £  =  local  width  of  trough 

6 .   Determination  of  Heat  Transfer  Rate 

Using  the  energy  balance  of  equation  (20)  with  the 
results  of  equations  (16)  and  (23),  the  heat  flow  for  a  fin 
element  of  width   d  x  and  length  H0   becomes 

dQ  =  U  dM  =  ^  ^  (ft^Ok+xs.n^as*  cos  gf  g^'25' 

9  ^  0^)4 


where  z*  =  z      5*  cos  o( 
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The  method  used  to  obtain  an  approximation  of  the  inside  wall 
temperature  gradient,   &  ,  was  to  equate  heat  flow  through 
the  fin  (equation  (25) )  to  that  which  could  flow  through  an 
element  having  width  EF  and  having  an  estimated  equivalent 
height  t*  (equal  to  t  plus  one-half  the  groove  height) . 

(26) 


(T»-T„-  ffi  _     (Z  k,  $f    Utjft  [(R^  k  s,n  <j>)  ^  cos  «  cos  $\x{?f 


t*i  3(^ 


This  resulted  in  a  transcendental  equation  for  6  which  was 
then  solved  using  an  iterative  scheme  at  each  incremental 
length  in  the  computerized  analytical  model  of  the  grooved 
condenser. 

The  mass  continuity  equation  for  any  x  location  of 
the  condenser  may  be  written  as 


1 


J^i     A  C    dM,       .  (27) 

Mtotai       \  *     dx    +  \  ~r^  d/ 

)0     dx  )o      dx 


where    £jji    is  obtained  from  equation  (25) ,  and 
dx 

dM 

L   is  obtained  from  equation  (24)  . 

dx 


The  first  step  in  evaluating  this  integral  is  to  divide  the 
condenser  into  a  large  number  of  incremental  elements  of 
length  A.X.   Equations  (24)  and  (25)  are  then  used  to  find 
the  mass  flow  rate  condensed  in  the  first   Ax  element.   This 
mass  flow  rate,  added  to  any  initial  flow  from  the  condenser 


29 


endpiece,  may  be  used  in  equation  (19)  to  solve  for  the 
trough  film  thickness,   &* ,  of  the  first  incremental  length. 
The  above  value  of    &*    is  then  used  in  equations  (24)  and 
(25)  to  find  the  mass  rate  condensed  in  the  second  incremental 
A  x  element.   The  above  process  is  repeated  until  the  entire 
condenser  length  has  been  traversed.   The  heat  transferred 
from  the  condenser  may  then  be  computed  for  each  incremental 
element  and  summed  over  the  length  using  equations  (24)  and 
(25)  or  obtained  from  the  product  of  total  condensate  mass 
flow  rate  and  the  latent  heat  of  vaporization  of  the  working 
fluid. 

E.   PARAMETRIC  STUDY  OF  INTERNALLY  GROOVED  CONDENSER 

The  theory  developed  for  the  internal  grooving  of  the 
condenser  section  was  intended  to  reduce  the  film  thermal 
resistance  and  increase  the  overall  heat  transfer  capability 
of  the  thermosyphon .   As  in  the  smooth  condenser  parametric 
analysis,  the  working  fluid  and  half  cone  angle  remained 
fixed.   For  the  grooved  geometry,  the  minimum  trough  width, 
£0 /  was  arbitrarily  chosen  to  have  the  same  dimension  as 
the  base  of  a  fin.   Plots  of  heat  transfer  rate  versus  fluid 
saturation  temperature  were  made  as  individual  values  of  the 

parameters  were  changed.   The  external  heat  transfer  co- 

2 

efficient  was  varied  from  1000  to  5000  Btu/hr  ft  °F  to  gen- 
erate one  family  of  curves.   Thermosyphon  speed  of  rotation 
was  varied  from  1000  to  3000  RPM  to  generate  another  family 
of  curves.   A  comparison  of  these  families  of  curves  was 
then  made  for  the  condenser  test  sections  involved,  stainless 
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steel  and  copper.   Also,  as  a  final  comparison,  external 
fins  of  the  type  discussed  in  the  smooth  condenser  study  were 
added  to  each  internally  grooved  test  section. 

F.   DISCUSSION  OF  ANALYTICAL  RESULTS 

In  actual  operation,  the  heat  transfer  rate  of  the  therm- 
osyphon  is  determined  by  its  application.   An  equilibrium 
saturation  temperature  at  which  the  thermosyphon  will  operate 
is  reached  for  all  other  operating  conditions.   Woodard's 
study  [6]  has  already  shown  water  to  be  the  most  satisfactory 
working  fluid. 

Figures   3-9   display  the  interrelationships  of  the 
previously  discussed  thermal  resistances  on  thermosyphon  per- 
formance.  While  an  increasing  rotational  speed  increases 
the  heat  transfer  for  any  fixed  external  h,  the  percentage 
change  becomes  larger  when  other  thermal  resistances  are 
lowered.   Proportionate  increases  from  1.25  to  1.70  are  ob- 
tainable upon  comparison  of  Figures   3   and   4  .   A  similar 
increase  for  the  internally  grooved  condenser  is  not  experi- 
enced in  Figures  6  and  7.   This  may  be  expected  since  off- 
setting effects  in  the  theory  make  heat  transfer  for  the 
grooved  condenser  essentially  independent  of  rotational 
speed. 

A  basic  explanation  is  that  a  higher  rotational  speed 
decreases  the  film  thickness  and  causes  an  increase  in  the 
heat  transfer  rate.   The  smaller  film  thickness,  however, 
increases  the  local  wall  temperature  and  causes  a  decrease 
in  the  heat  transfer  rate. 
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Similarly,  a  larger  external  heat  transfer  coefficient 
for  a  fixed  rotational  speed  is  seen  to  produce  proportionately 
larger  heat  transfer  rates  as  other  thermal  resistances  are 
reduced.   The  ratios  increase  from  1.05  to  7.0  for  the  com- 
parisons of  Figures  3,  5,  6,  and  8.   Figure  9  represents  the 
maximum  possible  increase  in  heat  transfer  under  the  given 
conditions  by  comparing  the  poorest  curve  of  Figure  3  to  the 
best  curve  of  Figure  6.   Of  note  is  the  fact  that  the  best 
curve  in( Figure  9  exceeds  the  entrainment  limit  of  the 
thermosyphon  as  discussed  by  Ballback  [3] .   Experimental  data 
of  Hewitt  et  al,  can  be  used  to  predict  the  entrainment  limit 
for  the  present  case  and  places  it  some  10-20%  below  the  con- 
denser limit  of  an  internally  grooved,  externally  finned, 
copper  test  section  [11] .   Since  the  entrainment  limit  is 
dependent  upon  the  square  of  the  condenser  diameter,  a  study 
could  be  performed  to  determine  that  diameter  at  which  the 
condenser  and  entrainment  limits  are  approximately  equal. 
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Figure  3   Heat  Transfer  Rate  vs.  Saturation  Temperature 
for  Internally  Smooth  Condenser 
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Figure  4   Heat  Transfer  Rate  vs.  Saturation  Temperature 
for  Internally  Smooth  Condenser 
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Figure  5   Heat  Transfer  Rate  vs.  Saturation .Temperature 
for  Internally  Smooth  Condenser 
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Figure  6   Heat  Transfer  Rate  vs.  Saturation  Temperature 
for  Internally  Grooved  Condenser 
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Figure  7   Heat  Transfer  Rate  vs.  Saturation  Temperature 
for  Internally  Grooved  Condenser 
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Figure  8   Heat  Transfer  Rate  vs.  Saturation  Temperature 
for  Internally  Grooved  Condenser 
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Fiqure  9   Heat  Transfer  Rate  vs.  Saturation  Temperature 
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III.   EXPERIMENTAL  PROGRAM 

A.   DESCRIPTION  OF  EQUIPMENT 

The  equipment  used  to  support  the  experimental  program 
was  the  same  equipment  used  by  Newton  [5]  and  Woodard  [6] . 
Some  changes  were  made  as  described  below  in  order  to  improve 
upon  data  accuracy  and  increase  overall  performance  of  the 
thermosyphon.   Figures  10  and  11  depict  the  experimental 
apparatus . 

1 .   Copper  Test  Section 

The  use  of  a  test  section  having  a  higher  thermal 
conductivity  was  previously  discussed.   Copper  was  chosen  as 
the  best  material  available  and  a  suitable  piece  of  pure 
oxygen-free  copper  was  purchased  for  construction  of  the  new 
thermosyphon.   A  new  condenser  and  evaporator  were  constructed 
The  condenser  was  machined  to  the  same  dimensions  as  the  1° 
stainless  steel  condenser.   Length  of  the  evaporator  had  to 
be  1/2  inch  shorter  than  its  stainless  steel  counterpart,  but 
all  other  dimensions  were  identical.   Stainless  steel  end- 
pieces  were  manufactured  for  the  evaporator  and  the  larger 
end  of  the  condenser  in  order  to  reduce  the  amount  of  axial 
heat  conduction.   The  evaporator  flanges  and  inconel  heater 
element  were  bonded  to  the  copper  section  in  a  hydrogen 
furnace  at  1850°F  using  a  50%  copper  -  50%  gold  brazing  foil. 
The  condenser  and  its  mating  flange  were  joined  using  high 
temperature  (1050°F)  silver  solder.   Assembly  of  the  copper 
evaporator  and  condenser  then  followed  the  guidelines 
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established  by  the  stainless  steel  thermosyphon.   New  bear- 
ings and  Garlock  seals  were  installed  on  the  new  test  section 
to  eliminate  some  of  the  noise  and  vibration  experienced  in 
rotation  of  the  stainless  steel  thermosyphon. 

2 .  Power  Supply 

Power  was  furnished  to  the  heater  coils  by  a  120  volt 
DC  motor  generator  having  an  automatically  regulated  voltage 
output.   The  use  of  parallel  connected  resistance  load  banks 
in  series  with  the  heater  element  gave  stable  heater  power 
levels  from  0.4  33  kw  to  a  maximum  of  8.0  kw.   To  meet  the 
higher  power  levels  that  an  improved  thermosyphon  design 
could  handle,  it  became  necessary  to  change  the  power  source. 
One  phase  of  a  three  phase  208  volt  AC  power  line  was  connected 
across  the  load  banks  and  heater  element.   Circuit  breakers 
were  set  to  trip  at  15  0  amperes     The  same  load  bank  arrange- 
ment was  used  to  regulate  the  heater  coil  power  level  through 
the  range  of  0.5  kw  to  20  kw.   Operation  of  the  thermosyphon 
with  the  new  power  source,  however,  created  unacceptable 
levels  of  random  noise  generation  in  all  eight  of  the  thermo- 
couples.  Brush  adjustments  and  attempts  to  ground  out  the 
noise  were  uniformly  unsuccessful.   Time  constraints  dictated 
a  return  to  the  DC  power  supply  which  continued  to  prove 
satisfactory  in  operation. 

3 .  Instrumentation  and  Coolant  Modifications 
a.   Thermocouple  Location 

To  obtain  a  better  temperature  profile  for  the 
condenser  wall,  six  thermocouples  were  used  in  an  uneven 
spacing  arrangement.   These  were  located  at  distances  of 
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0.0,  0.5,  1.25,  2.0,  3.0  and  5.3  inches  from,  the  flange  that 
mates  the  evaporator  and  condenser  sections  (Figure  10) . 
b.   No  Load  Heat  Generation 

Since  heat  was  being  generated  by  the  bearing 
and  Garlock  seal  and  transmitted  through  the  condenser  flange, 
it  became  evident  that  the  coolant  was  picking  up  some  portion 
of  this.   Both  copper  and  stainless  steel  test  sections  were 
run  under  no  load  conditions  to  determine  the  amount  of  heat 
generation  received  by  the  coolant. 

°l  coolant  =  mCp  AT  (28) 

where  AT  =  temperature  difference  of  coolant 
m  =  mass  rate  of  coolant  flow 

an  uncertainty  value  of  +  25%  was  assigned  to  this  data  due 
to  the  nature  of  the  source  and  the  measurement  techniques. 
Corrections  (as  shown  in  Table  2)  were  then  applied  to  the 
previously  obtained  data.   Even  though  the  flange  area  of  the 
copper  thermosyphon  was  insulated  with  low  conductivity 
silicone  rubber  during  assembly,  it  did  have  a  greater  amount 
of  no  load  heat  generation  than  the  stainless  steel  thermo- 
syphon.  It  is  felt  that  the  larger  number  of  thermocouple 
wires  strung  along  the  copper  condenser  imparted  this  in- 
crease in  energy  to  the  coolant. 
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RPM 

Q- STAINLESS 
*  (Btu/hr) 

Qz-COPPER 
(Btu/hr) 

700 

236+59 

^><c 

1000 

591+148 

1400 

394+97 

2100 

571+143 

866+216 

2800 

807+201 

1162+291 

Table  2   No  Load  Heat  Generation  for  Stainless  Steel  and 
Copper  Thermosyphons 

c.   Thermocouple  Attachment 

The  technique  for  welding  surface  temperature 
measuring  thermocouples  described  by  ref.  6  was  used  success- 
fully for  the  stainless  steel  condenser.   However,  this  pro- 
cedure did  not  work  for  the  copper  condenser  section.   Small 
grooves,  .010  inches  deep,  were  machined  in  the  wall  at  the 
desired  thermocouple  locations.   These  grooves  were  then 
filled  with  solder.   The  kapton  insulation  was  removed  from 
the  thermocouple  wire,  exposing  about  1/8  inch  of  the  invidual 
copper  and  constantan  wires.   The  bare  wires  were  then  cleaned 
and  the  tips  twisted  together.   This  junction  was  then  tinned 
with  the  same  solder.   A  small  oxygen-acetylene  torch  applied 
locally  to  the  inside  surface  of  the  condenser  would  melt  the 
solder  at  an  individual  thermocouple  location  and  the  tinned 
junction  could  easily  be  placed  at  the  desired  spot.   A 
strongly  bonded  thermocouple  junction  resulted  that  would  give 
the  average  surface  temperature  between  the  locations  that  the 
wires  penetrated  the  solder  surface.   The  thermocouple  wires 
were  then  tied  down  with  loops  of  copper  along  the  length  of 

the  condenser. 
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B.   EXPERIMENTAL  PROCEDURES 

The  experimental  procedures  involve  cleaning  of  the 
thermosyphon,  filling  with  the  working  fluid,  and  operating 
procedures. 

1 .   Cleaning  Procedures 

The  interior  surface  of  the  thermosyphon  had  to  be 
adequately  cleaned  to  insure  proper  wetting  by  the  working 
fluid.   If  the  stainless  steel  and  copper  test  sections  were 
not  cleaned  properly,  dropwise  rather  than  film  condensation 
would  be  obtained  with  water  as  the  working  fluid.   As  ex- 
plained in  Woodard ' s  study  [6],  the  theoretical  model  was 
based  on  film  condensation,  and  dropwise  condensation  results 
could  not  be  predicted.   Both  ethanol  and  freon  113  were 
found  to  have  excellent  wetting  characteristics  with  stainless 
steel  and  no  prescribed  cleaning  procedures  were  necessary. 
The  cleaning  procedure  described  in  ref.  6  was  satisfactory 
for  the  stainless  steel  thermosyphon  and  water.   Obtaining 
film  condensation  on  copper  required  different  cleaning 
procedures.   The  method  that  worked  was  a  modification  of  a 
procedure  used  in  preparing  copper  for  electroplating  [12] . 
The  following  is  the  appropriate  cleaning  procedure  for  the 
copper  thermosyphon: 

a.  Place  the  thermosyphon  in  the  vertical  position 
and  mount  the  fill  equipment. 

b.  With  end  windows  and  "0"  rings  removed,  fill  the 
thermosyphon  with  a  chlorinated  hydrocarbon 
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degreasing  compound  such  as  Trichlorethylene. 
Scrub  the  entire  surface  with  a  stiff  bottle 
brush.   (Avoid  prolonged  breathing  of  vapors) . 

c.  Drain  and  rinse  with  tap  water. 

d.  Fill  the  thermosyphon  with  a  mixture  of  sodium 
hydroxide  solution  and  200  ml  of  ethyl  alcohol. 
Warm  the  mixture  to  approximately  180°F.   Scrub 
with  a  stiff  bottle  brush. 

q.   Drain  and  rinse  with  tap  water. 

f.  Fill  the  thermosyphon  with  a  10-15  percent 
sulfuric  acid  solution  at  room  temperature. 
Scrub  with  a  stiff  bottle  brush.   Do  not  leave 
solution  in  thermosyphon  longer  than  5  minutes. 

g.  Drain  and  rinse  with  tap  water. 

h.   Rinse  three  times  with  distilled  water. 

i.   When  ready  to  fill,  prepare  the  end  windows  for 

installation,  drain  the  thermosyphon,  install  the 
end  windows  with  the  appropriate  "O"  ring,  gasket, 
and  bolt  ring.   Torque  the  bolts  to  30  lb-in. 
j .   Immediately  connect  the  fill  equipment  and  start 
the  vacuum  pump. 
2  .   Filling  Procedures 

Evacuation  of  the  thermosyphon  and  charging  with  the 
working  fluid  was  done  as  described  in  ref.  6.   Pumping  time 
for  the  vacuum  system  was  increased,  however,  to  allow  for 
complete  evaporation  of  any  liquid  remaining  in  the  thermo- 
syphon.  Longer  pumping  times  (10-20  hours)  and  use  of  the 
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diffusion  pump  gave  vacuums  of  less  than  100  microns 
pressure.   Care  was  taken  to  insure  the  mercury  manometer  on 
the  filling  apparatus  was  secured  after  the  initial  vacuum 
was  drawn.   The  mercury  vaporizes  easily  and  would  have 
contaminated  the  interior  surfaces  of  the  thermosyphon. 
3 .   Operational  Procedures 

Operating  procedures  remain  basically  the  same  as 
those  described  in  ref.  6.   Prior  to  applying  power  to  the 
heater  element,  the  thermosyphon  was  run  at  operational 
speeds  to  be  used  and  a  set  of  no  load  heat  generation  data 
taken  for  each  RPM. 

C.   EXPERIMENTAL  RESULTS 

1 .   Stainless  Steel  Condenser 

The  stainless  steel  thermosyphon  was  operated  quan- 
itatively  four  times.   De-ionized  distilled  water  was  used  as 
the  working  fluid  for  two  runs;  ethyl  alcohol,  and  freon  113 
were  used  for  the  other  two.   Limitations  on  operation  of  the 
thermosyphon  were  the  temperature  of  the  heating  coil 
(indicated  by  a  red  glow)  and  the  operating  pressure  in  the 
vapor  space.   Pressure  in  the  thermosyphon  was  not  allowed 
to  exceed  atmospheric  pressure  by  more  than  1-2  psi.   Such 
limitations  precluded  operation  at  low  RPM's  with  high  power 
levels.   As  indicated  by  Figures  12,  13,  14,  and  15,  water 
had  a  higher  operating  capacity  than  either  ethyl  alcohol  or 
freon  113. 

Visual  observations  were  made  through  the  end  window 
with  a  strobe  light.   In  all  operating  cases,  film  condensation 
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appeared  to  be  taking  place  since  the  condenser  surface  was 
mirror  smooth.   Only  in  the  water  run  of  3  May  1972  did  a 
spot  of  dropwise  condensation  appear  over  an  insignificantly 
small  portion  of  the  surface. 

2 .  Copper  Condenser 

Using  the  new  thermosyphon  of  copper/  a  data  run 
using  water  as  the  working  fluid  was  conducted  8  November  1972. 
Film  condensation  was  again  witnessed  over  all  but  a  very  small 
portion  of  the  condenser  surface. 

3 .  Coolant  Flov; 

During  the  data  run  with  the  copper  thermosyphon,  an 
attempt  was  made  at  modifying  the  external  heat  transfer  co- 
efficient.  Coolant  flow  was  increased  to  the  maximum  possible 
with  the  presently  installed  system.   Analysis  of  the  data 
showed  a  small  but  insignificant  increase  in  the  heat  transfer 
coefficient  from  1500  to  1650  Btu/hr-f t2-°F . 

A  different  modification  was  to  aim  one  of  the  coolant 
spray  nozzles  towards  the  evaporator  end  of  the  condenser 
(hottest  end) .   Data  analysis  again  showed  no  significant 
change  in  the  coefficient  or  heat  transfer  rate.   This  change 
of  coolant  pattern  did,  however,  flatten  the  condenser 
temperature  gradient  to  a  small  degree  as  shown  in  Figure  17. 

4 .  Discussion  of  Experimental  Results 

The  results  obtained  while  operating  both  thermosyphons 
are  plotted  in  Figures  12,  13,  14,  15,  and  16  together  with 
solutions  of  the  analytical  model  for  conditions  obtained 
during  operation.   The  two  water  runs  on  the  stainless  steel 
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condenser  establish  repeatability  since  the  same  cleaning 
and  filling  procedures  were  used  both  times.   Results  from 
these  two  runs  also  compare  very  closely  to  those  achieved 
by  Woodard  [6]  under  similar  conditions.   In  each  of  the  last 
two  cases,  a  better  thermosyphon  vacuum  than  Woodard ' s  runs 
was  achieved  prior  to  filling  (less  than  100  microns  Hg 
compared  to  4.5  mm  Hg) .   The  noted  difference  between  the 
experimental  results  and  the  analytical  curve  was  suggested 
to  be  a  result  of  the  presence  of  non-condensible  gases  in 
Woodard * s  study.   A  higher  degree  of  vacuum,  however,  would 
lessen  the  amount  of  non-condensibles  in  the  thermosyphons 
and  give  closer  experimental  and  analytical  results.   Such 
was  not  the  case  and  the  difference  cannot  be  readily  ex- 
plained.  It  may  be  speculated  that  a  minute  quantity  of 
non-condensible  gas  under  high  gravity  forces  would  tend  to 
have  the  same  blanketing  effect  as  a  greater  amount  of  non- 
condensible  gas. 

The  alcohol  and  freon  runs,  besides  having  a  lesser 
heat  transfer  capability  than  water,  had  the  unique  distinc- 
tion of  giving  experimental  data  that  exceeded  the  analytical 
solution.   This  trend  is  also  in  agreement  with  the  freon 
results  presented  by  Lee  and  Mital  [2] .   An  adequate  explana- 
tion cannot  be  deduced   other  than  to  note  that  constant 
magnitude  errors  for  all  runs  would  be  more  readily  apparent 
under  the  lower  heat  loads.   The  high  heat  transfer  capabil- 
ities of  water  could  overshadow  errors  that  become  significant 
for  alcohol  and  freon. 
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The  results  for  the  copper  thermosyphon  display  the 
same  trend  as  the  water-stainless  steel  runs.   Figure  16 
shows  that  the  heat  transfer  for  copper  was  less  than  stain- 
less steel  at  low  saturation  temperature  while  it  was  greater 
than  stainless  steel  at  higher  saturation  temperatures.   The 
results  for  each  still  do  not  differ  by  any  great  amount. 
This  was  expected  since  the  external  heat  transfer  coefficient 
and  internal  condensate  film  were  the  controlling  thermal 
resistances.   The  difference  between  experimental  and  theo- 
retical results  for  the  copper  test  section  is  much  larger 
than  the  differences  encountered  in  the  stainless  steel  runs. 
Analysis  of  the  data  has  indicated  that  the  thermocouple 
located  5.3  inches  from  the  evaporator  section  may  have  been 
reading  incorrectly.   The  temperature  distribution  over  the 
condenser  surface  was  generated  from  length  x=0  to  x=5.3 
inches  using  the  thermocouple  data  in  a  least  squares  poly- 
nomial and  then  from  x=5.3  to  x=9  inches  the  temperature 
distribution  was  given  a  constant  value  that  was  equal  to  the 
reading  of  the  thermocouple  located  at  5.3  inches.   Such  a 
temperature  distribution  is  illustrated  in  Figures  17,  18,  and 
19.   An  erratic  reading  for  the  thermocouple  located  at  x=5.3 
inches  would  therefore  give  erroneous  theoretical  results; 
even  accounting  for  the  curvature  at  the  start  of  the  theo- 
retical curve  for  the  copper  run  (Figure  16) .   Another  reason 
for  the  larger  difference  between  experiment  and  theory  may 
be  due  to  a  gradual  loss  of  thermosyphon  vacuum  (and  resultant 
introduction  of  non-condensible  gases) .   A  vacuum  of  100 
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microns  Hg  was  achieved  prior  to  filling  and  sealing,  and 
the  data  run  was  started  seven  hours  later.   Another  data 
run  on  the  succeeding  day  had  to  be  aborted,  however,  because 
the  maximum  capacity  of  the  thermosyphon  was  only  3  KW  of 
power.   Such  a  drop  in  heat  transfer  capacity  (from  8  KW  to 
3  KW)  with  time  can  only  be  explained  by  a  slow  loss  of  vacuum, 
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Figure  12   Thermosyphon  Heat  Transfer  Rate  vs.  Saturation 
Temperature,  3  May  1972  experimental  operation 
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Figure  13   Thermosyphon  Heat  Transfer  Rate  vs.  Saturation 
Temperature,  22  May  1972  experimental  operation 
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Figure  14   Thermosyphon  Heat  Transfer  Rate  vs.  Saturation 
Temperature,  2  June  1972  experimental  operation 
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Figure  15   Thermosyphon  Heat  Transfer  Rate  vs.  Saturation 
Temperature,  7  June  1972  experimental  operation 
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Figure  16   Thermosyphon  Heat  Transfer  Rate  vs.  Saturation 
Temperature,  8  November  197  2  experimental 
operation 
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Figure  18   Copper  Condenser  Outside  Wall  Temperature 
Profile  (Normal  Spray  Cooling) 
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Figure  19   Copper  Condenser  Outside  Wall  Temperature 
Profile  (Nozzle  inwards) 


60 


IV.   CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

1.  The  rotating  two-phase  thermosyphon  is  an  effective 
heat  transfer  device  capable  of  transferring  very  high  heat 
loads  at  high  RPM. 

2.  The  analysis  shows  that  performance  of  a  thermosyphon 
of  the  dimensions  used  in  this  study  is  strongly  influenced 
by  thermal  resistance  of  condensate  film  thickness  and  ex- 
ternal heat  transfer  coefficient. 

3.  The  performance  of  the  same  thermosyphon  is  influenced 
to  a  lesser  degree  by  the  thermal  resistance  presented  by  the 
wall  conductivity. 

4.  Effective  means  of  minimizing  the  three  thermal 
resistances  of  the  condenser  include  the  addition  of  external 
fins,  the  longitudinal  grooving  of  the  internal  surface,  and 
use  of  a  material  having  a  higher  thermal  conductivity. 

5.  The  performance  of  the  thermosyphon  may  be  degraded 
by  the  presence  of  non-condensible  gases  in  the  vapor  spaces. 

B.  RECOMMENDATIONS 

1.  Conduct  a  two-dimensional  conduction  analysis  in  the 
condenser  wall  to  determine  the  actual  temperature  distribution 
for  the  internally  grooved  condenser.   Implement  such  a 
solution  into  the  analytical  solution. 

2.  Investigate  heat  transfer  capabilities  for  internally 
grooved  condenser  having  other  than  triangular  profiles. 
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3.  Improve  the  external  cooling  method  by  addition  of 
fins  or  re-design  of  the  cooling  system. 

4 .  Conduct  a  systematic  study  of  the  effects  of  non- 
condensible  gases  on  thermosyphon  performance.   Devise  a 
method  to  determine  the  concentration  of  non-condensible 
gases  in  the  system. 

5.  Improve  the  condenser  surface  temperature  measurement 
of  the  copper  condenser  by  addition  of  more  thermocouples  and 
implementing  better  calibration  techniques  for  the  thermocouples 
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APPENDIX  A 


CALIBRATION 


1.   Calibration  of  the  Rotameter.   The  cooling  water  rota- 
meter was  calibrated  by  using  a  Toledo  scale  and  stopwatch 
to  measure  the  mass  rate  of  flow  of  coolant  flowing  into  a 
large  receptacle  placed  on  the  scale.   As  coolant  flowed 
into  the  container,  an  arbitrary  point  on  the  scale  was 
established  for  starting  the  stopwatch.   When  the  weight  of 
water  had  increased  by  40  pounds,  the  timer  was  stopped  and 
the  time  interval  recorded.   This  procedure  was  repeated 
three  times  at  each  flow  rate  and  the  average  time  interval 
was  calculated.   Data  was  taken  for  10,  20,  25,  30,  35,  and 
40%  full  rotameter  scale  flow  rates.   To  insure  that  the  flow 
duplicated  experimental  conditions,  the  coolant  was  allowed 
to  take  its  normal  flow  path  through  the  thermosyphon  cool- 
ing box.   The  drain  hose  from  the  cooling  box  then  emptied 
directly  into  the  weight  tank.   (The  scale  and  tank  were 
placed  in  the  gutter  to  maintain  proper  flow  in  the  drain 
hose.)   Data  was  plotted  (Figure  20  )  and  a  straight  line 
drawn  through  the  points.   The  equation  for  that  line  was 
found  to  agree  within  1%  of  the  flow  rate  measured  directly 
by  the  rotameter  scale  reading.   A  correction  factor  was 
calculated  and  the  actual  flow  rate  was  determined  from  the 
following  equation 

m  =  75  +  34  (%)  lbm/hr  (29) 

where  %  is  the  rotameter  setting. 
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2.   Calibration  of  Thermocouples.   Since  the  condenser  sur- 
face thermocouples  were  manufactured  in  place,  ordinary  means 
of  calibration  were  not  available.   Instead,  the  thermosyphon 
was  thermally  insulated  to  the  most  practical  degree  and  a 
source  of  steam  introduced  into  the  interior  of  the  condenser. 
After  the  condenser  had  warmed  up,  another  source  of  steam 
was  blown  directly  across  the  thermocouple  to  be  calibrated. 
In  this  manner,  it  was  determined  that  the  surface  thermocouples 
for  the  stainless  steel  thermosyphon  deviated  less  than  +  2.0°F 
from  the  steam  point.   A  larger  and  unacceptable  deviation 
occurred  for  the  thermocouples  on  the  copper  condenser.   It 
was  hypothesized  that  this  larger  difference  resulted  from  the 
high  thermal  conductivity  of  copper  and  inadequate  insulation. 
Calibration  for  the  copper  test  section  was  then  accomplished 
by  measuring  the  temperature  of  cooling  water  by  the  thermo- 
couples and  previously  calibrated  quartz  thermometers.   Agree- 
ment was  found  to  be  within  +  2.0°F. 


64 


1600  r 


1400  - 


1200 


1000 


05 

s 

PQ 


W 

O 
fa 


800 


600  - 


400  - 


200  - 


- 

0      / 

- 

0    / 

0    / 

O     t 

0     i 

• 

o    / 

/ 

• 

'                       » 

• 

• 

• 

10         20         30         40 

PERCENT  ROTAMETER  FLOW 


50 


Figure  20   Calibration  of  6.3  GPM  Rotameter 
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APPENDIX  B 
UNCERTAINTY  ANALYSIS 

Uncertainty  in  the  experimental  results  (shown  in 
Figure  12)  was  computed  in  accordance  with  the  procedure 
given  by  Holman  [13] .   The  heat  transfer  of  the  thermosyphon 
is  given  by 


Q  =  m  Cp  AT  -  qx 

where  Q    =  m  Cp  AT    (heat  from  no  load  generation) 

The  uncertainties  in  each  of  the  above  quantities  may  be 
expressed  as 


W.   =  +  17.0  lbm/hr 
m    — 


(20  to  1) 


Wr  =  +  00.001  Btu/lbm-°F   (20  to  1) 
<~p 

W.m  =  +  0.1       °C  (20  to  1) 

AT    — 

W  =  +   0.2  5   Q    Btu/hr   (20  to  1) 

The  uncertainty  in  heat  transfer  may  then  be  expressed  as 


W 


Q 


Q 


w. 

m 


m 


W 


+ 


AT 


W 


W 


+ 


+ 


AT 


C 


Using  the  3  May  1972,  stainless  steel  thermosyphon  run 
at  2800  RPM,  the  uncertainty  limits  become 
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Ts 

Q 

UNCERTAINTY 

94.0 

3431 

+  289 

6.8% 

112.9 

6374 

+  303 

4.2% 

130.0 

9180 

+  322 

O    .    /.  ~Q 

150.5 

12338 

+  348 

2.6% 

169.4 

15576 

+  380 

2.3% 

187.1 

18401 

+  411 

2.1% 

Uncertainty  in  no  load  heat  generation  proved  to  be  the 
dominating  factor.   As  shown  above,  its  effect  is  most 
noticeable  at  low  heat  transfer  rates. 
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APPENDIX  C 
TABULATED  DATA 

Water 
Stainless  Steel  Condenser 
Half  Cone  Angle  =  1° 


RPM 
Ts 

1000 

2000 

3000 

1 

H=1000 

100.0 

3994.5 

4558.8 

4845.9 

120.0 

6446.6 

7421.0 

7924.0 

140.0 

8877.8 

10270.8 

10994.0 

160.0 

11310.4 

13122.8 

14066.9 

180.0 

13741.9 

15972.0 

17137.0 

200.0 

16162.3 

18808.3 

20195.1 

H=1500 

100.0 

4532.9 

5303.6 

5714.5 

120.0 

7266.0 

8582.2 

9295.1 

140.0 

9969.6 

11841.2 

12861.4 

160.0 

12675.1 

15102.6 

16431.0 

180.0 

15381.0 

18361.4 

19997.6 

200.0 

18075.8 

21605.2 

23548.5 

H=5000 

100.0 

5539.9 

6807.5 

7548.8 

120.0 

8755.2 

10869.3 

12127.1 

140.0 

11923.7 

14890.5 

16669.3 

160.0 

15093.0 

18912.3 

21213.2 

180.0 

18266.6 

22934.3 

25755.2 

200.0 

21432.4 

26939.9 

30276.9 
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Water 
Copper  Condenser 
Half  Cone  Angle  =  1° 


RPM 

T 
s 

1000 

2000 

3000 

H=1000 

100.0 

4971.0 

5938.7 

6474.7 

120.0 

7917.3 

9550.4 

10470.0 

140.0 

10824.6 

13130.6 

14437.4 

160.0 

13731.3 

16708.2 

18401.6 

180.0 

16637.8 

20279.8 

22356.7 

200.0 

19531.7 

23870.8 

26287.7 

H=1500 

100.0 

5780.7 

7188.0 

8030.3 

120.0 

9099.8 

11429.8 

12846.4 

140.0 

12365.9 

15619.8 

17612.5 

160.0 

15630.4 

19805.8 

22373.7 

180.0 

18898.7 

23989.0 

27128.0 

200.0 

22158.7 

28152.9 

31856.2 

H=5000 

100.0 

7339.2 

9898.7 

11688.4 

120.0 

11279.7 

15316.0 

18162.5 

140.0 

15128.3 

20612.3 

24499.5 

160.0 

18967.4 

25890.7 

30814.7 

180.0 

22816.5 

31175.1 

37133.4 

200.0 

26669.7 

36456.8 

43443.2 
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Water 
Internally  Grooved  Stainless  Steel  Condenser 
Half  Cone  Angle  =  1° 

RPM=100  0 


H 
Ts 

100 
150 
200 

1000 

1500 

5000 

Fin  Half  Angle 
10° 

4253 
11124 
17991 

5295 
13811 
22322 

8251 
21451 
34616 

20° 

100 
150 
200 

4232 
11063 
17889 

5256 
13694 
22124 

8115 
21017 
33865 

30° 

100 
150 
200 

4226 
11040 
17848 

5232 
13617 
21989 

7990 
20628 
33202 

40° 

100 
150 
200 

4213 
10996 
17773 

RPM=3  0  00 

5194 
13502 
21797 

7827 
20159 
32425 

10°           100  4580  5763  9096 

150  12060  15125  23712 

200  19546  24495  38332 

20°           100  4564  5734  9007 

150  12014  15045  23455 

200  19473  24364  37902 

30°           100  4563  5723  8926 

150  12012  15009  23221 

200  19469  24305  37510 

40°           100  4557  5698  8802 

150  11994  14939  22871 

200  19438  24188  36930 
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Water 
Internally  Grooved  Copper  Condenser 
Half  Cone  Angle  =  1° 

RPM=1000 


H 

1000 

1500 

5000 

Fin 

Half  Angle 

Ts 

10° 

100 

5720 

7978 

20983 

150 

14914 

20764 

53680 

200 

24072 

33488 

85619 

20°  - 

100 

5692 

7899 

19699 

150 

14822 

20487 

49375 

200 

23910 

32996 

78218 

30° 

100 

5697 

7866 

18794 

150 

14815 

20342 

46893 

200 

23888 

32729 

74549 

40° 

100 

5709 

7843 

18142 

150 

14825 

20234 

45634 

200 

23894 
RPM=3000 

32532 

73243 

10° 

100 

6193 

8694 

22896 

150 

16266 

22698 

59634 

200 

26223 

36648 

96156 

20° 

100 

6180 

8660 

22346 

150 

16187 

22588 

57783 

200 

26157 

36461 

92912 

30° 

100 

6303 

8675 

21927 

150 

16240 

22608 

56342 

200 

26240 

36482 

90394 

40° 

100 

6234 

8701 

21530 

150 

16316 

22654 

55042 

200 

26360 

36544  • 

88172 
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TABULATED  EXPERIMENTAL  DATA 

8  November  1972 

Water 

Copper  Condenser 

Half  Cone  Angle  =  1° 

Standard  Cooling  Water  Nozzle  Arrangement 


RPM                    T            Q-EXPERIMENTAL  Q-THEORY 

1000              116.5  4323  8599 

151.9  9432  13428 

179.7  13186  17860 

2100              115.8  4520  12433 

141.0  9472  14048 

161.7  13284  18614 

214.3  24819  31344 

2800              115.5  5266  13442 

136.9  9747  14432 

155.4  13736  18744 

203.8  25094  31007 
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3  May  1972 
Water 
Stainless  Steel  Condenser 
Half  Cone  Angle  =  1° 
Standard  Cooling  Water  Nozzle  Arrangement 


RPM                T               Q-EXPERIMENTAL  Q-THEORY 

s 

700                95.1  3394  3364 

125.4  6199  7503 

150.5  9142  10661 
180.5  12321  14604 

207.0  15342  17660 

230.1  18128  20304 

1400               95.7  3295  4065 

117.5  6296  7468 

138.9  9298  10744 

164.0  12359  14404 

186.9  15734  17875 

206.0  18304  20695 

2100               95.3  3373  4301 

115.1  6355  7681 

132.4  9278  10405 

156.5  12300  14381 

174.6  15714  17533 
194.5  18303  20375 

2800               94.0  3431  4235 

112.9  6374  6576 

130.0  9180  10589 
150.5  12338  14196 
169.4  15576  17464 

187.1  18401  20108 
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22  May  1972 
Water 
Stainless  Steel  Condenser 
Half  Cone  Angle  =  1° 
Standard  Cooling  Water  Nozzle  Arrangement 


RPM 
700 


1400 


2100 


2800 


96.6 
123.0 
177.9 
203.5 

94.2 
116.7 
162.8 
184.0 
206.0 

94.0 
114.0 
155.1 
175.1 
194.5 


92 
115 

149 
166 
187 


Q-EXPERIMENTAL 

3448 

6309 
12344 
15185 

3447 

6327 
12401 
15282 
18671 

3485 

6463 
12420 
15281 
18651 

3602 

6423 
12458 
15358 
18650 


Q-THEORY 

3824 

7198 

14095 

17248 

3958 

7229 

14084 

17201 

20566 

4177 

7381 

14089 

17317 

20637 

4128 

7689 

13909 

16978 

20436 
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2  June  1972 

Ethyl  Alcohol 

Stainless  Steel  Condenser 

Half  Cone  Angle  =  1° 

Standard  Cooling  Water  Nozzle  Arrangement 

RPM  T  Q-EXPERIMENTAL       Q-THEORY 


700               88.0               925  954 

2185  1863 

4173  3558 

5118  4370 

984  1175 

2244  2063 

4232  3804 

5177  4646 

7303  6473 

1083  1282 

2264  2199 

4291  3928 

5275  4777 

7421  6748 

1122  1331 

2323  2214 

4331  3934 

5256  4717 

7460  6708 


Ts 

88. 

,0 

111. 

,2 

160. 

,9 

186. 

,2 

87. 

,5 

104. 

,6 

142. 

,4 

162, 

,0 

205, 

,9 

86. 

,6 

101, 

.3 

132, 

.7 

149, 

.3 

188, 

.2 

85, 

.3 

97, 

.9 

125 

.5 

138 

.9 

172 

.9 
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7  June  1972 

Freon  113 

Stainless  Steel  Condenser 

Half  Cone  Angle  =  1° 

Standard  Cooling  Water  Nozzle  Arrangement 

RPM  T  Q-EXPERIMENTAL       Q-THEORY 

700  82.7  807  470 

925  781 

1221  969 

1595  1278 

610  560 

984  929 

1259  1105 

1633  1439 

2086  1786 

630  602 

1023  984 

1279  1171 

1673  1537 

2145  1864 

610  600 

1043  1009 

1299  1180 

1693  1549 

2185  1898 


Ts 

82. 

,7 

96. 

,6 

105. 

,6 

121. 

,6 

81. 

,0 

93. 

.4 

99. 

,8 

112, 

,0 

125. 

,5 

80, 

,1 

90. 

,8 

96, 

,2 

107, 

,2 

117, 

.8 

78, 

.5 

88, 

.5 

92, 

.7 

102, 

.1 

111 

.8 
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